Numerical Visualization of Air Intake Induced by Free Surface Vortex 
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Free surface vortex control is vital in a pump sump system because the air absorbed by free surface vortex induc- 


es noise, vibration, and cavitation corrosion on the pumping system. In this study, the change of free surface vor- 


tex and air absorption in a pump intake has been investigated by the Volume of Fraction (VOF) method with 


steady multiphase flow model in order to represent the behavior of the free surface vortex exactly. The homoge- 


neous free surface model is used to apply interactions of air and water. The results show that air intake by the free 


surface vortex motion can be visualized using the iso-surface of air volume fraction. The vortices make an air 


column from the free surface to the pump intake. Also, it was found that the free surface vortex can be controlled 


by installing curtain walls. 
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Introduction 


Pumps are extensively used in cooling water system 
for electrical power plants. Good design of a pump sump 
is needed for pumping efficiency. However, sump design 
has been influenced by field conditions and economic 
costs. Poor sump structures lead to serious problems 
causing vortex. There are two types of vortex, free sur- 
face vortex in the approaching flow and submerged vor- 
tex induced on the walls. Free surface vortex is generated 
on the free surface and then absorbs air into the suction 
pipe of the pump. Such vortices induced by fluid flow on 
the free surface and walls can reduce pump efficiency or 
damage the impeller. 

Computational Fluid Dynamics (CFD) using numeri- 
cal approach is widely used to solve fluid flow and heat 
transfer. The results are worth and reliable, both in fluid 
flow and heat transfer analysis. Therefore many re- 
searchers have analyzed fluid flow near the pump intake. 
Nagahara et al. compared the PIV experiment of velocity 
near the submerged vortex and CFD results [1]. Rajen- 
dran et al. simulated three-dimensional turbulent flow 
including a vertical tailpipe inside a square channel and 
compared it with experimental results [2]. Padmanabhan 


and Hecker defined the shape of vortex which occurred 
near the tailpipe and classified it according to Froude 
number [3]. Okamura et al. used many CFD programs to 
solve the problem and analyzed the differences of the 
results [4, 5]. Chuang et al. calculated swirl angle from 
vorticities of cross section at the height of swirl meter [6]. 
Verhaart et al. constructed Anti Vortex Device (AVD) 
models and compared the velocity of inlet at suction pipe 
with respect to AVD dimensions [7]. Park et al. studied 
the change of free surface vortex according to various 
intake conditions [8]. 

In this study, the behavior of air intake caused by free 
surface vortex in the pump sump has been mainly inves- 
tigated using CFD. Numerical visualization technique 
with the aid of VOF method was adopted to express 
clearly the air absorption process induced by surface 
vortex. Also, the function of a curtain wall to control the 
free surface vortex is studied. 


Numerical analysis 


Governing equation 
The governing equations used to analyze the flow field 


of pump sump are the continuity equation and the Navi- 
er-Stokes equation with standard k-e turbulence model. 
The equations are as follows: 
Continuity equation: 
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where u; (i=1, 2, 3) is the velocity component in the i 
direction; f; is the quality forces in the i direction; p is 
pressure, p is fluid density, u is fluid viscosity;,G,, the 
production of k, is given by 


G, =u, [Side (5) 
Ox; Ox; } Ox; 
Where, uis the eddy viscosity: 
k2 
= pC, F (6) 


The k-e model coefficients are adopted as following: 

C,=0.09, o=1.0, 051.3, Ca=1.44, Cy=1.92, 

In this study, the homogeneous model for multiphase 
flow analysis is applied. The density and physical quan- 
tity are represented by the volumetric ratio of multiphase 
flow on each of the working fluids. 


p, U, Fare defined as follows 
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Each fluid y1, y2, ..., N, is the total number of phases. 
Volume fraction, density, velocity, diffusion coefficient of 
each fluid are Ya, Pas Uas Ta respectively. a represents the 
total phase. 


Computational domain and boundary conditions 


Fig.l shows the computational domain of pump sump 
model used in the numerical analysis. This study is refe- 
renced by the model proposed by Okamura et al. [5]. 
Three-dimensional modeling, grid generation and com- 
putational model calculations were performed on this 
model. All dimensions were based on the outside diame- 
ter (D) at inlet of suction pipe. 

Analysis software used for this purpose is ANSYS 
CFX-16.2, the number of grids used are 0.8~1.1 mil- 
lion nodes, after checking the grid independency, 
steady calculations were performed. Homogeneous 
multiphase flow model is applied in order to imple- 
ment the multi-phase flow of air and water, expressed 
as the volumetric ratio of the volume of fluid. On the 
inlet boundary condition, the pressure condition p=pgh 
was applied - using CFX user function (expression) - 
to express the pressure gradient according to water 
level. Outlet boundary is applied for the flow rate 
condition. 

Flow rate and water level are closely related to the 
vortex creation. If the flow rate is too high, then the wa- 
ter level should be increased in order to prevent vortex 
occurrence. The relationship between flow rate and water 
level is expressed as 

S=D(1.0+2.3F p) (10) 
where S is the minimum water level, Fp is Froude num- 
ber. Minimum water level S is the standard height to 
prevent the free surface vortex without any additional 
constructions [9]. 

The parameters for the water levels and flow rates are 
given in Table 1. For case 1, free surface vortices occur 
usually [5]. Hence, case 1 was adopted to describe the 
numerical visualization of air intake induced by free sur- 
face vortex. In case 2, curtain wall is used to control the 
free surface vortex. The dimensions and location of the 
curtain wall are based on the HI Standard recommenda- 
tions. Case 2 is used to check the effects of the curtain 
wall. 

Table 1 Operating conditions of numerical model 


Flow rate (m/h) Water level (mm) Curtain wall 


Case 1 135 3.75 D No 


Case 2 135 3.75 D Yes 


Results and discussion 


Free surface vortex visualization by the Volume Frac- 
tion method 

For case 1 and 2, the real fluid flow is 135 m°/h, but 
Fp is based on 60 m°*/h. Therefore, cases 1 and 2 are both 
severe water level conditions and hence it is easy to 
create free-surface vortices. 
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Fig. 1 Computational domain of pump sump model for nu- 
merical analysis 


Fig. 2 shows the streamline near suction pipe for case 
1. The Free-surface vortex is identified at both sides of 
the suction pipe (Fig 2(a)) and it can be seen that the flow 
near the suction pipe inlet and free surface are connected 
to each other (Fig 2(b)). This flow is called free-surface 
vortex and it consists of mostly water and small parts of 
air bubble. For visualization of air intake induced by 
free-surface vortex, air volume fraction has been 
checked. 

Fig. 3 indicates the vorticities and iso-surface of air vo- 
lume fraction from 5% to 0.1%. When the volume fraction 
ratio of air is more than 5%, vortex is not visualized; when 
the volume fraction is at 1%, the free surface vortex is 
observed and if less than 1%, it can be clearly observed 
from the free surface towards the inlet pipe of the pump 
making the air-water column as shown in Fig 3 (c) and (d). 

In other words, when using the homogeneous model, 
the free surface vortices are expressed at the region of 
specific volume fraction by changing the air volume 
fraction factor. Free surface vorticities are observed for 
the case when air or water volumetric ratio is 5-0.1% air 
and 95-99.9% water. 

Fig. 4(b) shows the side view streamline around the 
suction pipe with curtain wall for case 2. Flow near the 
free surface is blocked by the curtain wall and it is sepa- 
rated into upstream and downstream unlike Fig. 4(a). The 
downstream flow flows to inlet of the suction pipe di- 
rectly and rapidly. Due to the curtain wall, upstream flow 
makes recirculation zones under the free surface. These 
two flows are independent and the direction of movement 
is opposite to each other. Therefore, the vertical flow 
creating column from free surface to the inlet of suction 
pipe does not exist. The result can be shown through 
vortex visualization by air volume fraction. 
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(b) Side view 


Fig.2 Streamline near suction pipe for case 
Effect of curtain wall for air intake 


Curtain wall is a kind of AVD that is installed to con- 
trol free-surface vortices in severe cases. Curtain walls 
are of various shapes, but the vertical type is widely used. 
The curtain wall used in case 2 is also of the vertical type 
and its dimension specifications are determined from the 
HI standard [9]. However, the dimension of 1.0D is used 
as the depth of the curtain wall, instead of 0.5 of the 
recommended dimension from HI standard, as Pyo Y.S. 
proposed that 1.0D is a more reasonable depth in order to 
control the free surface vortex [10]. 

Fig. 4(b) shows the side view streamline around the 
suction pipe with curtain wall for case 2. Flow near the 
free surface is blocked by the curtain wall and it is sepa- 
rated into up and downstream unlike Fig. 4(a). The 
downstream flows to inlet of the suction pipe directly and 
rapidly. Due to the curtain wall, upstream flow makes 
recirculation zones under the free surface. These two 
flows are independent and the movement of direction is 
opposite to each other. Therefore, the vertical flow creat- 
ing column from free surface to the inlet of suction pipe 
does not exist. The result can be shown through vortex 
visualization by air volume fraction in Fig. 5 (a). 

Fig. 5 shows the effect of the presence of the curtain 
wall for air intake with 0.1% air volume fraction com- 
paring vorticities and iso-surfaces. For the case of no 
curtain wall as shown in Fig. 5 (a), it can be seen that, the 
vortex flow move towards the inlet pipe of the pump 
from the free surface, thus creating air water column. In 
the case of Fig. 5 (b), there is no strong vertical flow 
creating air-water column from free surface to the inlet of 
the suction pipe due to the presence of a curtain wall. 
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(a) Air volume fraction 5% 
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(c) Air volume fraction 0.5% 
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(b) Air volume fraction 1% 
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(d) Air volume fraction 0.1% 


Fig.3 Air intake induced by free surface vortex 
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(a) No curtain wall 


Fig. 4 Comparison of streamline between no curtain wall and with curtain wall 
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(b) With curtain wall 
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(a) No curtain wall 
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(b) With curtain wall 


Fig.5 Effect of curtain wall for air intake 


f i A numerical visualization technique using iso-surface of 
ps air volume fraction has been applied for a better under- 

— standing of the air intake mechanism. The numerical re- 
¢ sults show that using the iso-surface of air volume fraction, 

the air intake by the free surface vortex motion can be vi- 

sualized. The vortices make an air column from the free 

surface to the pump intake. Also, it was found that the free 

surface vortex can be controlled by installing curtain walls. 


Conclusions 


The locations and behavior of free surface vortices 
have been studied in a pump sump using multi-phase 
flow simulation with steady condition. Numerical visua- 
lization technique with the aid of VOF method was 
adopted to express clearly the air absorption process in- 
duced by surface vortex. The conclusions are as follows: 

1) Air intake induced by free surface vortex can be 
clearly visualized by expressing the streamlines and the 


iso- surface of constant air volume fraction. 


2) It reveals that free surface vortex can be controlled 


by installing the curtain wall. 


3) In order to remove free surface vortex, the vertical 
flow created from free surface to the inlet of suction pipe 
has to be controlled. Curtain wall installation is a good 
choice for removing this vertical flow. 
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